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ABSTRACT 

We present observations of the molecular gas in the nuclear environment of 
three prototypical low luminosity AGN (LLAGN), based on VLT/SINFONI AO- 
assisted integral-field spectroscopy of H 2 1-0 S(l) emission at angular resolutions 
of ~ 0.17". On scales of 50-150 pc the spatial distribution and kinematics of 
the molecular gas are consistent with a rotating thin disk, where the ratio of 
rotation (V) to dispersion (<r) exceeds unity. However, in the central 50 pc, the 
observations reveal a geometrically and optically thick structure of molecular gas 
iy/a < 1 and Nh > 10 23 cm -2 ) that is likely to be associated with the outer 
extent of any smaller scale obscuring structure. In contrast to Seyfert galaxies, 
the molecular gas in LLAGN has a V/a < 1 over an area that is ~ 9 times 
smaller and column densities that are in average ~ 3 times smaller. We interpret 



^ased on observations at the European Southern Observatory VLT (082.B-0709 and 084.B-0568). 
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these results as evidence for a gradual disappearance of the nuclear obscuring 
structure. While a disk wind may not be able to maintain a thick rotating struc- 
ture at these luminosities, inflow of material into the nuclear region could provide 
sufficient energy to sustain it. In this context, LLAGN may represent the final 
phase of accretion in current theories of torus evolution. While the inflow rate is 
considerable during the Seyfert phase, it is slowly decreasing, and the collisional 
disk is gradually transitioning to become geometrically thin. Furthermore, the 
nuclear region of these LLAGN is dominated by intermediate-age/old stellar pop- 
ulations (with little or no on-going star formation), consistent with a late stage 
of evolution. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: kinematics and 
dynamics — infrared: galaxies 

1. Introduction 

For three decades, the "Unified Model" of Active Galactic Nuclei (AGN) has dominated 
discussion. Its most popular (and extreme) version assumes that AGN Types are distin- 
guished only by viewing angle, with the broad lines in Type-2s blocked from our line-of-sight 
(LOS) by a highly-inclined molecular /dusty torus. This premise is supposed to hold for 
all luminosities, from low luminosity AGN (LLAGN) to powerful quasars (Antonucci 1993). 
More realistic formulations of the torus model involve additional variables such as covering 
factor or AGN luminosity (see the discussion on Elitzur 2012). Direct evidence for this 
model comes from spatially-resolved interferometric observations of thermal dust in a few 
prototypical AGN (e.g., Jaffe et al. 2004; Tristram et al. 2007). However, these observations 
do not provide sufficient information on the actual geometry (e.g., position-angle coverage 
necessary to constrain shape/orientation) and physical properties (e.g., dynamics) of the 
obscuring medium. 

By means of Adaptive Optics (AO) integral-field spectroscopy, Hicks et al. (2009, here- 
after H09) showed that the molecular gas in the central tens of parsecs of Seyfert galaxies 
relates directly to the largest structures associated with the obscuring torus, as predicted by 
clumpy torus models (e.g., Nenkova et al. 2002, 2008; Schartmann et al. 2008): it is in a 
rotating disk-like distribution, has a high velocity dispersion relative to rotation (V/a < 1), 
and is optically thick (see also Miiller- Sanchez et al. 2006, 2009). HCN measurements of 
Seyfert galaxies suggest that the nuclear dense gas also has a large dispersion (Sani et al. 
2012). However, it remains unclear how this vertical structure can be supported. Some 
models suggest that the torus may be an evolving structure driven by accretion processes 
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operating on scales of tens of parsecs such as nuclear starbursts or disk- instabilities (e.g., 
Vollmer et al. 2008; Wada et al. 2009; Hopkins et al. 2012). When the inflow rate is low, the 
thickness and opacity of the torus decrease. Alternatively, b y considering both th e broad-line 



region and the torus to be different aspects of a disk-wind, lElitzur &: Hoi (120091 ) argue that 
the torus thickness depends on AGN accretion. They predict that below a given luminosity 
threshold, a wind with a significant column density (Nh) can no longer be sustained. In 
either case, the geometrically and optically thick structure of molecular gas in LLAGN is 
expected to disappear. 

Observationally, there is an on-going debate on the existence of an obscuring structure in 
LLAGN0. The non-detection of an IR-bump in the SEDs of LLAGN (IFernandez-Ontiveros et al 



20121 ) provides evidence for unobscured nuclei, as do the compact variable central UV sources 
observed in both types of LLAGN (e.g., Maoz et al. 2005). Furthermore, X-ray studies of 
nearby galaxies have found a correlation between Nh and X-ray luminosity (e.g., Zhang et 
al. 2009). However, some LLAGN do exhibit high Nh (e.g., NGC 3169; Terashima & Wil- 
son 2003). The detection of polarized broad- lines and ionization bicones in some LLAGN 
provides additional support for the presence of an obscuring torus (Barth et al. 1999; Pogge 
et al. 2000). 

In an effort to characterize for the first time the properties of the molecular gas in the 
central regions of LLAGN, we have observed three prototypical LINERs/LLAGN (Table 1) 
using AO-assisted i^-band integral-field spectroscopy. The 2.12 /zm H 2 1-0S(1) emission line 
is used to analyze the two-dimensional (2D) distribution and kinematics of the gas within 
the central 150 pc. These techniques have already been applied to study the properties of H 2 
in NGC 1097 (H09; Davies et al. 2009). Here we bring those data together with new data 
on three additional objects (§2 and §3) to probe the putative obscuring medium in LLAGN. 



2. Observations and Data Processing 



The galaxies analyzed here were observed with AO and SINFONI at the VLT in two 
separate runs between March 2009 and March 2 010. SINFONI provide s spectra for every 
pixel in a contiguous 2D field with 64 x 32 pixels (jEisenhauer et al.ll2003l ). The pixel scale in 
all observations was 0.05" x 0.1", resulting in a field-of-view of 3.2" x 3.2". The spectral range 



1 The bulk of the LLAGN population reside in low-ionization nuclear emission-line regions (LINERs; Ho 
2008). Since the bolometric luminosity threshold for separating low- and high-luminosity AGN is roughly 
defined as a few 10 42 erg s _1 (Ho 2009), some LLAGN reside in Seyfert galaxies. Strictly speaking, we study 
LINERs/LLAGN here, but refer to them just as LLAGN. 
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was similar for all runs (1.95-2.45 /xm) with typical spectral resolution R~4000 (FWHM~ 70 
km s -1 ). For each galaxy, the nucleus was used as NGS by the AO-module. The exposure 
times for NGC 1052, NGC 2911, and NGC 3169 were 180, 20, and 60 minutes, respectively. 
The FWHM angular resolution, measured from the spatially unresolved non-stellar contin- 
uum in if— band, is ~ 0.17" in NGC 1052 and NGC 3169, corresponding to 15 and 16 pc, 
respectively, and ~ 0.23" (45 pc) in NGC 2911. Data reduct ion was performed using the 



SINFONI data reduction package SPRED (lAbuter et al.ll2005l ). Nearby standard stars (A 



or B-type) were observed close in time to the science frames and used for telluric correction 
and flux calibration. In addition, flux calibration was cross-checked with VLT/NACO and 
2MASS photometric data in 1" — 4" apertures. 

The resulting spectra show several emission lines of H 2 and stellar absorption features 
(Fig. 1). The 2D properties of H 2 1-0 S(l) were extracted using the code LINEFIT (Davies 
et al. 2009). This procedure corrects automatically for instrumental broadening. Typical 
uncertainties in the kinematic measurements are between 5 — 20 km s _1 . 



3. Results and Analysis 

The SINFONI data reveal the morphology and kinematics of the molecular gas in un- 
precedented detail (Fig. 2; see H09 for the maps of NGC 1097). Similar morphologies are 
observed in the three sample galaxies. The flux distributions show bright central structures 
with elongations from Northeast to Southwest and morphologies suggesting filamentary for- 
mations. While in NGC 3169 (morphological type S), the photometric position-angle (PA) 
of H 2 is similar to the PA of the stellar continuum, in NGC 1052 and NGC 2911 (type E4 
and SO, respectively), the PA of H 2 is strongly misaligned (by ~ 70 — 80°) with that of 
the stars and the galactic major axis, indicative of an external origin of th e gas, probably 



being acquired through merging, as has been proposed for elliptical galaxies (IMorganti et al. 



20061 ). 



The velocity fields of H 2 exhibit azimuthal symmetry, a zero-velocity axis (or kinematic 
minor axis) almost perpendicular to the kinematic major axis (which is also the photometric 
major axis in all cases), and small twists of the iso- velocity contours. Taken together, these 
features provide strong qualitative evidence for rotation in a disk, but with the presence of 
radial flows in the nuclear region. In NGC 1052 and NGC 2911, the gas has a kinematic 
major axis that is not aligned with the stellar rotation axis (Fig. 2), providing strong support 
to the hypothesis of an external origin of the molecular gas. 



While the H 2 velocity fields indicate rotating disks, the dispersion maps suggest that 



-5 - 



these disks are fairly thin in the outer r > 0.4" (a ~ 20 — 30 km s _1 ), but remarkably 
thick in the vicinity of the AGN (a increases up to ~ 150 km s _1 ). In all cases, a also 
increases in a direction almost perpendicular to the H 2 kinematic major axis (in NGC 3169 
this is observed only in the Southeast, Fig. 2). Our experience with previous integral-field 
observations of other galaxies, reveals that increases in dispersion in particular directions 
surroundi ng the AGN, are usually related to the presence of a passing radio-jet and/or 
outflows (iMiiller-Sanchez et al.ll201ll ). Radio and X-ray observations of N GC 1052 and NGC 



3169 reveal extended emission in the direction o: 



maps of NGC 2911 flWrobel et al. 



1984 



Hummel et al.l 1 19871: Terashima fc Wilson 2003). A similar situation is observed in radio 



the enhancements in a ( Kadler et al.ll2004 



Condon et al.lll99ll ). Thus, the regions of enhanced 



a- values (delineated by dashed lines in the dispersion maps, Fig. 2), are probably due to 
radio-jets or AGN-driven outflows. 

For quantitative modeling we adopted an axisymmetric rotating disk generated by the 
code DYSMAL (as in NGC 1097, Davies et al. 2009; 2011). For determining the disk 
parameters and their uncertainties, we adopted an exponential mass distribution of outer 
radius 1.6". We then varied inclination, PA, total mass, black hole (BH) mass and radial 
scale length (Ri/ e ) to minimize x 2 m the residual velocity and dispersion maps. The regions 
perpendicular to the kinematic major axis containing high a-values (the area inside the 
dashed lines in the dispersion maps, Fig. 2), which are likely of non-gravitational nature, 
were not included in the fit. The smearing due to the instrumental resolution is taken into 
account by convolving the inclined model with the Gaussian point-spread function of the 
appropriate width. Input errors in the x 2 evaluation were the la-fit uncertainties of the line 
profiles. All uncertainties given for the derived parameters are 68% confidence. 

The best-fitting models are presented in Table 1 and illustrated with one-dimensional 
cuts along the kinematic major axis (Fig. 3). The models provide a good fit to the observed 
LOS-velocity and a curves, but the relatively large nuclear dispersion (at r < 0.4") cannot 
be accounted for by the models. The dominant kinematic component in all cases is well 
described by circular motion. However, the velocity residuals maps of NGC 1052 and NGC 
2911 show non-circular velocities of up to ~ 60 km s -1 within 1.5" of the nucleus. Since 
the purpose of this investigation is to derive the best-fitting velocity field and study the 
bulk properties of H 2 , we do not discuss the residuals here (including the regions where 
the outflows are observed). A detailed analysis of the kinematics in these two galaxies will 
be presented in another publication. Overall, the observed properties of the nuclear H 2 in 
these LLAGNs indicate a disk-like distribution with significant rotation and high velocity 
dispersion in the central r < 0.4". 



-6 - 



4. Discussion 

On scales of r > 50 pc (~ 100 pc in NGC 2911) the properties of the molecular gas are 
consistent with a thin disk model (Fig. 3). Inside this radius, the dispersion increases faster 
than the models, and within r < 25 pc (~ 50 pc in NGC 2911 but we use 25 pc for comparison 
purposes) the ratio V/a is < 1, implying that the gas on these scales must be supported 
by rotation and random motions within a spheroidal structure (e.g., a geometrically thick 
disk). In addition, the gas in this region is optically thick. From the estimated dynamical 
masses, and assuming a typical gas fraction (f g ) in the centers of galaxies between 1 — 
10%§ (H09; Hopkins et al. 2012), we obtain an average N H of 1.7 x 10 23 cm" 2 (Table 
1). In NGC 1052 a nd NGC 3169 the estimated Nfj is consistent with that derived from 



X-ray spectroscopy (ITerashima &: Wilson! 120031 ; iKadler et al.l 120041 ) . Taken together, these 



characteristics suggest a rotating, geometrically and optically thick disk, and are consistent 
with those expected for the molecular torus in unification models. A similar conclusion was 
reached for the nuclear H 2 in NGC 1097 (H09). Therefore, we refer to this structure as "large- 
scale torus", keeping in mind that the obscuring material probed by mid-IR interferometry 
of Seyferts has a size of a few pc or less (e.g., Jaffe et al. 2004), and that this is likely 



respo nsible for other functions such as collimation of ionization cones (iMuller-Sanchez et al. 



20111 ). Nevertheless, the contribution of the large-scale torus to AGN obscuration may be 



considerable, as suggested by current clumpy torus models (e.g., Schartmann et al. 2008). 

The large-scale torus in LLAGN appears to be different from the one in Seyfert galaxies. 
This conclusion is apparent from the ratio of azimuthally-averaged V/ a as function of radius 
(Fig. 4). For comparison purposes, we also plotted the average of this function for six Seyfert 
and two starburst galaxies observed with SINFONI on the same spatial scales (r ~ 150 pc). 
On average, for LLAGN, V/a is low (0.8 ± 0.3) at r = 25 pc and high (4.6 ± 1.2) at 100 
pc. Thus random motions are significant within r < 25 pc, despite the rotation of the gas. 
In contrast, starburst galaxies exhibit a high V/a at any radius consistent with a thin disk 
(an almost uniform dispersion map with values < 30 km s _1 , i.e. no significant thickening 
at the center), and Seyfert galaxies show low V/a at radii up to 100 pc, suggesting that the 
geometrically thick disk in these objects is more extended than in LLAGN. The obscuration 
in LLAGN is probably also lower than in Seyferts. The average 1-0 S(l) luminosity of the 
LLAGN in our sample within r = 25 pc is ~ 1.5 x 10 4 L Q , approximately 3.2 times smaller 



2 Th e gas mass can also be est imated from the luminosity of 1-0 S(l) using an appropriate conversion 
factor ( Miiller- Sanchez et al. 20061) . We used the same factor as H09 of 430 M Q /L Q (which gives f g ~ 10% 
for Seyferts), and derived an average f g ~ 3.2% in the central 50 pc of our sample of LLAGN. This is the 
value used throughout the paper. CO measurements of LLAGN also suggest f g of a few percent (Raluy et 
al. 1995; Tan et al. 2008). 
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than the average found for Seyfert galaxies (H09). If the H 2 luminosity reflects gas mass 
(rather than differences in excitation), then the LLAGN in our sample have gas masses (and 
in consequence Nh) on average ~ 3.2 times smaller than Seyferts (H09). This is consistent 
with CO measurements of local AGN (including NGC 1097 and NGC 1052; Raluy et al. 
1998; H09), which indicate that there is more gas in the centers of Seyferts than in LLAGN. 

Fig. 4 clearly shows that the molecular gas in LLAGN is in an intermediate state between 
an extended large-scale torus (a thick disk) and no large-scale torus at all (a thin disk). We 
interpret this structure as a large-scale torus in the process of gradual disappearance. This 
conclusion is supported by the theory of torus evolution proposed by Vollmer et al. (2008). In 
this model, a torus passes through 3 different phases dictated by an external mass accretion 
rate. Initiated by a massive gas infall, a turbulent thin disk is formed in which starbursts 
occur. Such disks, characterized by the production of several star clusters and high star 
formation rates (SFR > 10 M yr" 1 ), are observed in local starburst galaxies (Muller- 
Sanchez et al 2010, Fig. 4). Once supernovae remove the intercloud medium, a collisional 
thick disk of dense gas clouds is formed with a decreasing, but still high accretion rate. The 
collisional thick disk also forms stars, but with an efficiency of ~ 10% when compared to the 
previous stage (SFR ~ 1 M yr _1 ). The Seyferts in the SINFONI subsample (Fig. 4) are 
interpreted to be in this second phase (H09). The mass inflow rate is slowly decreasing, and 
the collisional disk is gradually transitioning to become geometrically thin. The SFR reaches 
a minimum, and the nuclear region is dominated by intermediate- age/old stellar populations. 
Based on the observed properties, the LLAGN in our sample are in this transitional phase 
(Figs. 1 and 4). When the inflow rate has significantly decreased, the collisional torus either 
completely disappears (there is very little gas left) or becomes thin and transparent as the 
circumnuclear disk in the Galactic Center. 

The proposed scenario reveals a possible Starburst — » Seyfert — » LINER/HII — > 
LINER — >■ Passive Galaxies evolutionary sequence of AGN activity and BH growth. The 
initial phase of activity of a galactic nucleus is characterized by vigorous star formation. 
After ~ 10 Myr, significant accretion onto the BH triggers the Seyfert phase. This highly 
active phase may be modulated by feeding and feedback processes. The final stages of 
accretion are observed in LINERs, in which massive BHs have already formed. This picture 
is supported by recent observational studies of local galactic nuclei (e.g., Constantin et al. 
2009; Schawinski et al. 2010), and numerical simulations of the origin and evolution of 
luminous AGN (e.g., DiMatteo et al. 2005; Hopkins et al. 2006). 

The kinematic signatures of outflows detected in the dispersion maps suggest that a 
disk-wind may still be present. In the wind scenario, the main channel for release of energy 
may be switching at low luminosities from AGN-driven outflow to radio-jets. Thus, the 
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regions in which we detect high a could be interpreted as gas from the ISM shocked by a 
radio-jet. However, it seems unlikely that this mechanism is also responsible for the high 
dispersion measured in the central 50 pc (at least not in its most part). First, the circular 
symmetry of these regions (extending out to r ~ 25 — 50 pc) is not consistent with that 
expected from a collimated jet. In addition, in some objects, there is no clear evidence of a 
radio-jet (e.g., NGC 1097; Orienti & Prieto 2010), and yet a low V/a is observed (Fig. 4). 
Although a contribution from the disk-wind to the large-scale torus seems unlikely, it could 
significantly contribute to sustain any smaller scale torus (the threshold for a Mbh = 1 x 10 8 
M is L bol > 2 x 10 40 erg s -1 ; Elitzur & Ho 2009). The dust sublimation radius in LLAGN 
with L bol = 10 42 erg s" 1 is r sub = 0.4(L bo i/10 45 ) 1 / 2 pc ~ 0.012 pc (Nenkova et al. 2008). 
The outer radius of the "small-scale" torus is then R = 50 — 100r sn b = 0.6 — 1.2 pc (Nenkova 
et al. 2002; Tristram & Schartmann 2011), which, at our resolution, is unresolved. These 
results suggest that the torus comprises at least two components on different scales each of 
which fulfills different roles. On scales of ~ 15 pc and larger, our results reveal a diffuse 
and clumpy envelope which represents the reservoir of gas necessary to feed the AGN, and 
possibly contributes to nuclear obscuration. At smaller scales, IR interferometry provides 
evidence for a compact torus which is likely responsible for most of the obscuration and 
collimation. 

Our integral-field data confirm the presence of an obscuring structure of molecular gas 
in the central 50 parsecs of LLAGN, which indicates that in these objects there is potentially 
enough material to fuel the AGN. However, it appears to be disappearing as indicated by 
the lower column densities and sizes compared with those of Seyfert galaxies. Particularly 
interesting would be a study of the gas kinematics in AGN having lower luminosities and/or 
Eddington ratios than the ones presented here. In these objects either there is no detectable 
nuclear warm-H 2 structure (as in M87 and some other LLAGN observed with AO; Gebhardt 
et al. 2011, A. Barth priva te communication) , or it is thin and transparent as predicted by 



the torus evolution model ( IVollmer et al.l 120081 ) . 
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Fig. 1. — X— band nuclear spectra of LLAGN extracted from our SINFONI data in r = 0.3" 
apertures. The prominent emission lines are from H 2 (marked at the top of the figure), 
the strongest being 1-0 S(l) at 2.12 /xm. Several stellar absorption features (like the 12 CO 
bandheads or Nal), characteristic of evolved stellar populations, are also noticeable. All 
K— band ionization lines (including Br7) are absent in these objects. NGC 1097 exhibits a 
similar spectrum, but Br7 is marginally detected (Davies et al. 2007). Since Br7 is normally 
considered a tracer of young stars, the marginally/non-detection of this line indicates little 
(as the very young and compact starburst in NGC 1097, Davies et al. 2007) or no ongoing 
star formation. 
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Fig. 2. — Two-dimensional maps of H2 1-0 S(l) flux distribution, LOS-velocity, and velocity 
dispersion (from left to right), for each galaxy observed with SINFONI. Contours delineate 
the if— band continuum emission. The position of the AGN (the peak of non-stellar con- 
tinuum at 2.2 (Urn) is marked with a cross. The dashed line in the middle column shows 
the orientation of the stellar kinematic major axis measured from the SINFONI data. In 
all cases this is consistent with the major axis of the galaxy. The two dashed lines in the 
dispersion maps delineate approximately the regions where the radio and/or X-ray emission 
is observed in each galaxy. Rejected pixels in the velocity and dispersion maps are those 
with a flux density lower than 5% of the peak of H 2 emission. 
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Fig. 3. — Velocity and dispersion distributions of the molecular gas along the H 2 kinematic 
major axis of the sample galaxies (see H09 for the curves of NGC 1097). Asterisks represent 
LOS- velocity measurements and open diamonds dispersion values. Continuous curves denote 
the best-fitting exponential disk model to the velocity datapoints, and dashed curves the 
same models but for a. In all panels vertical error bars represent la uncertainties in the 
measurements. 
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Fig. 4. — Azimuthally-averaged V/a versus radius for each galaxy as indicated by the legend. 
The dotted horizontal line shows the point where V/a = 1. The typical error of the V/a 
estimates is 35%. The average V/a as a function of radius was obtained for Seyferts from 
the d ata of NGC 3783 and NG C 7469 (H09), NGC 2992 (Friedrich et al. 2010), NGC 
6814 jMuller-Sanchez et alJl^OllI ), NGC 1566 and NGC 3081 (this work); and for starburst 
galaxies from NGC 253 and NGC 4945 (this work). 
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Table 1: Summary of basic data for LLAGNs and quantities derived from the H 2 1-0 S(l) kinematics. 



"Spectroscopic classification and Distance (D) according to Ho et al. (1997), except for NGC 1097 (H09). 
L2: LINER without a broad Ha line, LI. 9: LINER with weak broad Ha emission, LI: LINER with broad 
components of both Ha and H/3 lines. L2s in the sample show convincing evidence of be ing LLAGN. The two 



have a high brightness-temperature radio core, which is also a hard X-ray point source (jTerashima fe Wilson 
20031: rWrobel et al.lll984t lEllis et al.ll2006l) . 



h Bolometric luminosity estimated using the X-ray luminosity as Lboi = 15L2-iokeV (Ho 2009), except for 
NGC 1097, which was estimated from the SED analysis (Prieto et al. 2010). X-ray luminosities taken from 
Ellis et al. (2006) and Terashima & Wilson (2003). 

c Black Hole mass derived from the dynamical models, except for NGC 1097 (Davies et al. 2007). 
^Dynamical mass estimated taking into account the velocity dispersion of the gas as Md yn = (V 2 + 3a 2 )r/G, 
where r is the radius from the nucleus and G is the gravitational constant. 

e Hydrogen column density derived assuming the typical 1 — 10% gas mass fraction (this interval corresponds 
to the lower and upper limits). Quoted values refer to 3.2% gas fraction (see text for details). 

^Dynamical mass obtained from the kinematic modeling, except for NGC 1097 (H09). This is consistent with 
the estimate derived taking into account the velocity dispersion at this radius within the uncertainties. 



Note. — Quoted uncertainties are 1 standard deviation of the fit, except for Nn, which is delimited by 
the assumed gas fraction interval. 



